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ABSTRACT 

The astrophysical site of the r-process is still uncertain, and a full exploration 
of the systematics of this process in terms of its dependence on nuclear proper- 
ties from stability to the neutron drip-line within realistic stellar environments 
has still to be undertaken. Sufficiently high neutron to seed ratios can only 
be obtained either in very neutron-rich low-entropy environments or moderately 
neutron-rich high-entropy environments, related to neutron star mergers (or jets 
of neutron star matter) and the high-entropy wind of core-collapse supernova 
explosions. As chemical evolution models seem to disfavor neutron star mergers, 
we focus here on high-entropy environments characterized by entropy S, electron 
abundance Y e and expansion velocity V exp . We investigate the termination point 
of charged-particle reactions, and we define a maximum entropy Sf ina i for a given 
V exp and Y e , beyond which the seed production of heavy elements fails due to the 
very small matter density. We then investigate whether an r-process subsequent 
to the charged-particle freeze-out can in principle be understood on the basis of 
the classical approach, which assumes a chemical equilibrium between neutron 
captures and photodisintegrations, possibly followed by a /3-flow equilibrium. In 
particular, we illustrate how long such a chemical equilibrium approximation 
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holds, how the freeze-out from such conditions affects the abundance pattern, 
and which role the late capture of neutrons originating from /3-delayed neutron 
emission can play. Furthermore, we analyze the impact of nuclear properties 
from different theoretical mass models on the final abundances after these late 
freeze-out phases and /3-decays back to stability. As only a superposition of astro- 
physical conditions can provide a good fit to the solar r-abundances, the question 
remains how such superpositions are attained, resulting in the apparently robust 
r-process pattern observed in low metallicity stars. 

Subject headings: abundances, charged-particle reactions, high-entropy wind, nu- 
cleosynthesis, r-process, supernova 



INTRODUCTION 



A rapid neutron-capture process (r-process) in an explosive scenario is traditionally be- 
liev ed to be responsible for the nucleosynt hesis of about half of the heavy elements above 



Fe ( jBurbidge et al.l 119571 ; ICameron 1 119571 ) . While in recent years the high-entropy wind 
(HEW) of core-collapse supernovae has been considered to be one of the most promising 
sites, hydrodynamical simulations still encounter difficulties to reproduce the astrophysical 
conditions under which this process occurs. Therefore, a model- independent approach, i.e., 
the s o-called 'Waiting- Point" approximat i on fW P), has been utilized for many years (see, 



e.g., ( ICowan et al.lll99ll ; Kratz et al.lll993l . 120071 )). This rather simple model, with the basic 



assumptions of an Fe-group seed, an (n, 7) — (7, n)-equilibrium for constant neutron densi- 
ties n n at a chosen temperature T over a process duration r and an instantaneous freezeout, 
has helped to gain steadily improved insight into the systematics of an r-process in terms 
of its dependence on nuclear-physics input and astrophysical conditions. Taking a specific 
seed nucleus (for convenience often 56 Fe was employed), the solar r-process pattern with its 
three pronounced abundance peaks at A=80, 130 and 195 can be reproduced for a variation 
of neutron number densities n n and a given temperature T. Whether the solar r-process 
residuals N r0 ~ N — N s<3 is fully reproduced in each astrophysical event, i.e., whether each 
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In realistic astrophysical environments with time variations in n n and T, it has to be investi- 
gated whether at all and for which time duration r the supposed (n, 7) — (7, n)-equilibrium of 
the classical approach will hold and how freeze-out effects change this behavior. In general, 
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late neutron captures may alter the final abundance distribution. In this case, neutron cap- 
ture reactions will be important. Also /3-delayed neutrons (/3dn) can play a role in forming 
and displacing the abundance peaks after freeze-out. 

An example of a more realistic astrophysical environment is the HEW expelled from newly 
formed (hot) neutron stars in core-collapse s upernovae, wh i ch continue to release n eutrinos 
after the supernova shock wave is launched ( jMeyer Ill993l ; iQian fc Woosley Ill996l ). These 
neutrinos interact with matter of the outermost proto-neutron star layers which are heated 
and ejected in a continuous wind. T he late neutrino flux also leads to moderately neutron- 
rich matter ( IQian &: Woosley 1119961 ) via interactions with neutrons and protons and causes 
matter ejection with high entropies. However, from the beginning problems were encoun- 
tered to at tain entropies sufficiently high in order to produce the heaviest r - process nuclei 
(see, e.g., (Takahashi et al.lll994t IWoosley et al.lll994j ; iThompson et al .11200 It IWanajo et al. 
200ll : iTerasawa et alJbooS )). Recent hydrodynamic simulations for core-collapse supernovae 



support the idea that these entropy constraints may be fulfilled in the late phase (after the 
initial explosion) when a reverse shoc k is forming (IFryer et al.l 120061 ; I Ar cones et al.l 120071 ; 



Burrows et all 120071 : Ijanka et al.l 120071 ) . 



The question is whether such high entropies o ccur at times with sufficiently high tem- 
peratures when an r-process is still underway (IWanajo I l2007t Kuroda et al.l 120081 1 . Ex- 
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)). In the 



present calculations we focus on (a) the ef fects of different nu clear physics input (mass models 
FRDM (Finite Range Droplett M odel. jMoller et al.l Il995h . ETFSI-1 (Extended Thomas- 



Fermi with Strutinsky Integral), (lAboussir et al.l 119951 ). a versio n with quench i ng of shell 



closures proportional to the distance from stability ETFSI-Q. dPearson et al.l 119961 ). the 
mass formula of Duflo & Zuker (DUFLO-ZU KER). (iDuflo fe Zuker I 119961 ) and a recent 



Hartree-Fock-Bogoliubov approach (HFB-17), ( Goriely et al. 2009^ ). and (b) a detailed un- 



derstanding of the r-process matter flow far off stability, in particular testing equilibria, 
freeze-out and effect of delayed neutron capture. To investigate these effects we have applied 
a full network containing up to 6500 nuclei and the corresponding nuclear masses, cross 
sections and /3-decay properties. 

Starting in the 1990's, in various publications discussing r-process nucleosynthesis, neutrino 
interactions were predicted to significantly affect the final r-process abundances. Neutrinos 
can affect r-process environments in four different ways at different stages: (a) by determin- 
ing the neutron-richne ss of matter (Y P ) y ia neutrino or antineutrino captures on neutrons 
and protons (see, e.g., (jMeyer et al.lll992l )). (b) by destroying a-particles by neutral current 
interactions leading to higher seed production, i.e., reducing the neutron to seed ratio (see, 
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e.g., (IMeyer Ill995l )). (c) by speedin g up the flow to heavy nuclei by charged curren t reactions 
mimicing fast /3~-decays (see, e.g., (IFuller fc Meyer 111995c iQian fc Woosley Ill996l ) or finally 
(d) by acting via neutral-current (spallation) reactions to possibly "filling up" unde rabun 
danc e s in the low-mass wings of the A~130 and neutrino- induced fission (see, e.g., (iQian 



20021 ; iKolbe et al.l 120041 ) . ( b) acts like starting an r- process with a higher Y e , (c) was found 
to be rather unimportant (IFreiburghaus et al.lll999l), in (d) the effect s of /3-delayed neutron 



emission seem to be dominant ( jKratz et al. 



2001 



Pfeiffer et al.ll2001al ). and ^-delayed fission 



was also found to be unimportant. Thus, our calculations are based on trajectories for the 
matter density p(t) and the temperature T(t), originating from the HEW expansions. An 
extended parameter study of the r-process has been performed in terms of entropy S, elec- 
tron abundance Y e and expansion velocity V exp of the hot bubble, the latter quantity being 
related to the expansion timescale r exp . 



2. THE COMPUTATIONAL FRAMEWORK: HIGH-ENTROPY 
EXPANSIONS AND NUCLEAR PROPERTIES 

In the absence of self-consistent hydrodynamical models which also result in r-process 
conditions permitting to produce matter all the way up to the third peak and U and Th, we 
continue to utilize parametrized calculations in order to survey the dependence on nuclear 
properties and highlight a detailed understanding of the nucleosynthesis processing which has 
so far not been fully analyzed. As mentioned in the introduction, many of such calculations 
have been performed in the past which displayed final abundance distributions. We intend 
here to focus on monitoring for which time intervals equilibria are obtained, for which nuclear 
mass ranges this applies, and what are the effects of the emission and recapture of /3-delayed 
neutrons on the (final) abundances. 

We fo llow the description of adia batically expanding homogeneous mass zones as already 



utilized in IFreiburghaus et al.l ( 119991 ). Different mass zones have different initial entropies so 
that the overall explosion represents a superposition of entropies. The electron abundance 
Y e , the entropy S and the expansion velocity V exp are parameters representing the degrees 
of freedom related to the expansion properties for the thermodynamic evolution and the 
nucleosynthesis in an adiabatic expansion. 
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2.1. Thermodynamics 



If the pressure per unit volume by relativistic particles (photons, electrons, positrons ..) 
is in excess of the total pressure due to non-relativistic particles, the entropy is dominated 
by radiation. This occurs at high temperatures and moderate to small matter densities. As 
high temperatures also ensure a nuclear statistical equilibrium composition, it is possible to 
choose an arbitary, but sufficiently high, initial temperature, e.g., T 9 ~ 9 (T 9 = T/10 9 K). If 
the thermal energy fc^T is larger than the energy of the rest mass of electrons, we have also 
to consider the contibutions by electrons and positrons. The radiation dominated entropy 
is given by 5 7 = (4/3)aT 3 /p, where 5 7 is the entropy per unit mass and a is the radiation 
constant. Making use of relations for ultrarelativistic fermions, the electron and positron 
contributions are given by S e + >e - = {7 /4)S 1 = (7/3)aT 3 / p. The total entropy per unit mass 
is therefore 



S Siy ~\~ *Sg+ ) g— ~ Ob 



11 T 3 
-a- 

3 p 



(1) 



which is a factor (11/4) larger than S T After the temperature has decreased to values where 
the corresponding energies become comparable or smaller to the rest mass of electrons, the 
contributions of electrons and positrons can be neglected. The electrons become nonrela- 
tivistic and the positrons cease to exist. To consider both e xtreme situations (p ure radiation 
or radiation plus ultra- relativistic electrons and positrons), IWitti et al.l (119941 ) introduced a 
very close approximation for the entropy 

7 m) 



S — 5*-v 



i + 



with 



f(T 9 ) 



4* 



T 2 
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(2) 



(3) 



The fit function f(T$) varies between and 1. If we write the entropy S in units of kg per 
baryon and the density p in units of 10 5 g cm -3 (p 5 ), the entropy 5* can be expressed as 



S 



1.21^ 



T 3 
Ps 



(4) 



Since the expansion of the hot bubble proceeds adiabatically (TV 1 / 3 = const), the time 
evolution of the temperature is given by 

R 



T 9 (t) = T 9 (t = 0) 



Rq + Vexp t 



(5) 



where R is the inital radius of an expanding sphere and V exp is the expansion velocity of 
that sphere. From Eq. H] one can deduce the matter density p 



PM = 1-21 



T 3 

1 9 
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(6) 
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Making use of a typical initial value Ro = 130 km and a sufficiently high initial temperature 
Tg = 9, which ensures nuclear statistical equilibrium to consist only of neutrons and protons, 
we have performed calculations for a large grid of expansion velocities (1875, 3750, 7500, 
15000 and 30000 km/s, for a given S) and a large grid of entropies (5 < S < 415). It 
should be noted that for the lowest entropies of this grid (S < 15) the above treatment is 
probably not valid, as the radiation dominance is not fulfilled. However, for the purpose of 
an r-process survey this parametrization seems still useful. 

Fig. [TJ shows the time evolution of T 9 (t) and p 5 (t), both normalized to 1, for a selected 
expansion velocity of V exp = 7500 km/s. These normalized quantities are identical for 
different entropies as long as V exp is kept constant. 

The electron abundance Y e = YliZiYi (with proton number Zi and abundance Yi = Xi/Ai 
[mass fraction over mass number]) is equal to the averaged electron or proton to nucleon 
ratio < Z/A > and thus is a measure of the proton to neutron ratio in the HEW. This is of 
key importance for the r-process after the freeze-out of the charged-particle reactions. A 
neutron- rich wind means that Y e < 0.5. For an inititial Tg = 9 only free protons and neutrons 
exist, leading to 1 = X n + X p = Y n + Y p = Y n + Y e . Therefore, a given Y e defines the initial 
proton and neutron abundances Y p = Y e and Y n = 1 — Y e . The initial value of Y e depends 
ultimately on all weak interactions (including neutrino interactions) with the available free 



in the early p] 


rase of ( 


(Frohlich et al. 


2006a 


b 



are expected to experience Y e < 0.5 and enable the onset of an r-process. In order to test 
such conditions we have chosen a large grid of Y e (from 0.40 up to 0.499) for our calculations. 



2.2. Nuclear Networks and Nuclear Properties 

2.2.1. The Charged-P article Network 



The nucleosynthesis calculations during the early (hot phase in the) expansion of 
the wind until char ged-particle freeze - out, were performed with the Basel nucleosynthe - 
sis code (see, e.g., ([Thielemann et al.l Il996l ; iBrachwitz et al.l l2000t iFrohlich et al.l l2006al ); 
however, not including neutrino interactions) making use of the nuclear input as described 
there. Some new experimental rates were added since, but for the majority of two-particle 
rates and their inverse reactions - especially those for unstable nuclei - the predictions 



of iRauscher fc Thielemann I (120001 ) (rate set FRDM) were used. The code includes all neu- 
tron and charged-particle induced reactions as well as their inverse reactions, /3-decays, elec- 
tron captures, /3dn emission and particle transfer reactions. The charged-particle nuclear 
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network considered for all simulations of this paper in the early expansion phase, before the 
onset of the r-process, is shown in Table [TJ 



2.2.2. The r-Process Network 



The r-process calculations after charmed-particle freeze-out made use of the updated 



version of the code as documented in iFreiburghaus et al.l (119991 ). which includes neutron- 
induced reactions and /3-decay properties as well as a simple formulation for fission. The 
neutron capture and /3-decay rates were up dated and modified in order to study the im- 
pact of a variety of mass models (F RDM ( Moller et al. 1995h . ETFSI- 1 ( Aboussir et al. 



1995h . ETFSI-Q dPearson et al1ll996h. and HFB-17 rtGorielv et aZIboogh ) and a mass for- 
mula DUFLO-ZUKER ( iDuflo fc Zuker Ill996l ). The size of the network was chosen specific 
to each mass model because the neutron dripline is located at different mass number A 
for each mass model. Specific care was taken to include (3dn emission and their re-capture 
also during r-process freeze-out. Therefore, temperature-dependent neutron capture rates 
had t o be introduced. For these purp oses both the rate sets based on FRDM and ETFSI- 



Q of iRauscher fc Thielemann I (l2000f ) were used for the experimentally unknown rates of 



neutron-rich nucle i. Details on the nuclear spectr oscopy used in the calculation of these 
rates are given in IRauscher fc Thielemann I (j200lh . To study the impact of further mass 
predictions by ETFSI-1, DUFLO-ZUKER and HFB-17, the effective neutron separation en- 
ergies S n used in the computation of the photodisintegration rate, were modified according 
to the mass predicitions. 

The method described here, making use of a splitting between the charged-particle network 
(until charged-particle freeze-out) and a (fast) r-process network for the neutron capture 
phase thereafter, is valid as long as the initial phase (which includes charged particle reac- 
tions) is not producing nuclei beyond the limit of the charged particle network (Z=46). This 
is the case for all conditions in our parameter s urvey of Y K , S , and expansion timescale. Cal- 
culations with very short expansion timescales (IMeyer 1120021 ) differ, but we did not consider 
them here, as it has not been established that such conditions reproduce the abundance 
pattern of heavy r-process nuclei, which is characterized by /3-decay properties acting on 
longer timescales. 



3. EXPLOSIVE BURNING AND CHARGED-PARTICLE FREEZEOUT 



At initially high temperatures, 9~ 9, the baryonic matter mainly consists of free neu- 
trons and protons. Then follows the recombination of the free nucleons into a-particles. 
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The formation of a-particles proceeds until the temperature drops to values around 9~ 6 
when the a-particles dominate the nuclear abundances. The maximum abundance of the 
a-particles which can be attained at that time, depends on Y e . If Y e < 0.5 and all (initial) 
protons are consumed to form a-particles, their abundance will be 

Y Y 

Y = — = — (7) 

1 max 2 2' 

and the corresponding mass fraction is 

X m ax ^Y e . (8) 

The mass fraction and the corresponding abundance of the remaining free neutrons is 

X n = Y n = 1 - 2Y e . (9) 

Table [2] shows the maximum mass fraction X max and the corresponding X n before any seed 
can be synthesized, and Fig. [2] shows the behavior of the a-particles for a selected Y e - value 
and six different entropies. 

Between 9~ 6 and 3, depending on the matter density, a complete or incomplete recombi- 
nation of the a-particles can take place. Thereby, the still existing neutrons will help (a) to 
overcome the bottle-neck regions at A=5 and 8 and (b) to build neutron-rich seed nuclei. 



3.1. The a-Rich Freezeout 

When during the expansion of the hot bubble the temperature drops below 9~ 6, the 
a-particles start to react with each other to form heavier nuclei via 12 C, either following the 
triple-a -reaction 

3a ->• 12 C, (10) 
or, if neutrons are available, the reaction 

a + a + n^ 9 Be, (11) 

followed by 9 Be(a,n) 12 C. 

Both Eqs. [10] & [11] depend on the square of the matter density, and Eq.fTTl depends in addition 
on Y e . They represent bottle-necks for the subsequent nucleosynthesis towards heavier nuclei 
and are the reason for the well-known difference between a normal and an a-rich freeze-out of 
charged-particle reactions. Further reactions of possible importance like 3 He+ 4 He — > 7 Be+7, 
7 Be + 4 He -> n C + 7, 4 He + 2 H ^ 6 Li + 7, 6 Li + 4 He ^ 9 Be + 1 H, and 9 Be + 4 He 
C + 7 are already included. In addition, test calculations with Fynbo's triple-alpha reaction 
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rates (IFynbo et al.ll2005l ) have been performed for two entropies S = 100 and 250 kg/baryon, 
although in several recent publications the data analysis with respect to the "missing" 2 + 
rotational-band memb er of the 12 C Hoyl e state and its astrophysical consequences have been 
questioned (see, e.g., ( iFreer et al.l 120091 )). We also have tested the NACRE triple-a rates 
for the above two entropies, and no significant differences in the seed contributions were 
observed at the charged-particle freeze-out. For high matter densities, i.e., low entropies 
(S oc T 3 /p), the three-body reactions of Eqs. [10] and [11] will remain effective, leading to 
a total recombination of a-particles and free neutrons. The final abundances at charged- 
particle freeze-out close to 9~ 3 are dominated by iron group elements. For a moderate 
neutron excess or 0.4 < Y e < 0.5, the free neutrons will be consumed completely and no 
subsequent rapid neutron capture (r-process) will b e possible. This is referred to as a normal 



freeze -out of charged-particle reactions (see, e.g., (jWoosley et al.lll994l ; iFreiburghaus et al. 



19991 ) and references therein). 



In the case of lower matter densities, i.e., higher entropies, the reactions of Eqs. [TO] and [TT1 
cease to be effective for further recombination of a-particles towards heavier nuclei. The 
degree to which such a production of heavier nuclei is prohibited is a gradual function of S 
and Y e , ranging from a-particle mass fractions of a few percent to close to 100 X max = 200 Y e 
% (e.g., for y e =0.45 the maximum amount of the a-particles is 90%) and accordingly small 
amounts of heavy "seed" nuclei (see Fig. EJ). 

In an a-rich freeze-out the full nuclear statistical equilibriu m (NSE) among all nuclei breaks 
down into quasi-equilibrium (QSE) subgroups (see e.g., (jffix fc Thielemann I fl999h ). For 
a strong a-rich freeze-out and moderately neutron-rich Y e , the dominant abundances of 
heavy (seed) nuclei during the charged-particle freeze-out can be shifted to mass numbers 
between 80 and 110, thus overcoming the shell closure at N=50. The ratio of neutrons 
to heavy seed nuclei at this point is a function of Y e (measuring the neutron excess) and 
entropy (determining the ratio of heavy nuclei to a-particles and indirectly also the amount 
of neutrons consumed in these heavy nuclei). This can be seen from Fig. El and will be further 
elaborted in the following subsections. The neutron to seed ratio is a measure of the amount 
of neutron captures on seed nuclei and the resulting average mass number after an r-process, 
and indicates therefore the strength of the r-process. The fact that the seed abundances are 
dominated by nuclei with mass numbers between 80 and 110, is one reason that an r-process 
in the HEW can be faster than in classical calculations, which typically starts below N=50. 



3.2. Freezeout Timescale 



We define the expansion time T exp between our arbitrarly chosen initial temperature 
T 9 = 9 and its decrease to T 9 = 3 as the timescale for the charged-particle reactions and 
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as the time when the freeze-out of such reactions occurs within the expanding and cooling 
wind. According to Eq. [5] one gets 

_ 2^o n9 N 

T exp — y j K 1 ^) 
'exp 

where Rq is the initial radius of the hot bubble and V exp is its expansion velocity during the 
successful supernova explosion. For instance an expansion velocity of 7500 km/s corresponds 
to an expansion timescale of 35 ms. 



3.3. Seed Distribution 

After an a-rich freeze-out, matter which surpassed the bottle-neck above A=4 is pre- 
dominantly accumulated in nuclei of the Fe-group or beyond. At sufficiently high entropies 
these nuclei are shifted to the mass numbers in the range 80 < A < 110. If free neutrons are 
still available, neutron captures on those seed nuclei can proceed after the freeze-out of the 
charged-particle reactions, thus leading to a subsequent r-process. For a given Y e , (V exp , S)- 
combinations can be chosen which result in a fixed neutron to seed ratio (Y n /Y see( j) and a 
unique distribution of seed nuclei at freeze-out (independent of the expansion timescale for 
such choices). This is shown in Fig. H] for a selected Y e - value with a constant neutron to 
seed ratio and for 4 different expansion velocities. This means that for an r-process with 
a given strength for a specific Y e , i.e., with a given initial neutron to seed ratio, both the 
initial temperature (T 9 w 3) and the seed distribution are always the same. However, in 
order to get the same r-process ejecta, i.e., in some sense a universal r-rprocess, similar ex- 
pansion velocities (timescales) are required. The reason for this is that the r-process path for 
a given strength also depends on the temperature profile which is relate d to the expansion 



fr 

timescale. Similar to other authors (lArcones et al.l 120071 ; iWanajo 1 120071 ). we use the terms 



"hot" or "cold" r-process dependending on the average temperature experienced when the 
neutron captures occur (which enables photodisintegrations to be effective or not). For in- 
stance, a small expansion velocity (V exp < 3000 km/s) leads to a "hot" r-process which is 
completely governed by an (n, 7) — (7, n)-equilibrium. On the other hand, higher expansion 
velocities lead to a "hybrid" r-process, which is first equilibrium dominated at high temper- 
atures (Tg > 1), then followed by a non-equilibrium phase at lower temperatures (T 9 < 1), 
i.e., the "cold" r-process where neutron captures and /3-decays compete. In Fig. [5] we show 
the dependence of the seed composition for the same wind characteristics (i.e., the same 
strength and timescale), but for different electron abundances Y e . In contrast to the quite 
robust seed distribution of Fig. HJ here we observe for a constant V exp a significant variation 
with Y e , with the major differences at the highest seed masses above A« 90. Since the mass 
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fraction of a-particles increases with decreasing Y ej the amount of the ejected heavy material 
(X heavy = 1 — X a ) is different for the same r-process strength, provided that the expelled 
matter contains roughly the same amount of mass per entropy interval. 



3.4. Termination of the Charged-Particle Reactions 

There is a natural upper limit of the r-process- relevant entropies, i.e., the last entropy 
which still yields a significant amount of seed nuclei. We label those maximum entropies as 
S final- Beyond Sfi na i the recombination of the a-particles, i.e., Eqs. [10] & HH fail completely 
and the yields at such entropies only consist of a-particles, free neutrons which decay to 
protons, and traces of 2 H , 3 H, 3 He, 12 C and 16 0. This is similar to the big bang nucleosyn- 
thesis where the dominant protons and the less abundant neutrons recombine to roughly 
75% protons, 24% 4 He and traces of 6 Li and 7 Be. Table |3] shows an example of the yield 
of such a boundary entropy beyond Sf ina i. In the limiting case of disappearing protons, the 
abundance and mass fraction of a-particles and the neutrons are given by Eqs. [7] to IHJ and 
are shown for some values of Y e in Table |2j 

Table [4] finally shows the results of simulations for two expansion velocities and a selection 
of Y" e -values. The second and the sixth column contain the initial entropies (Si n iu a i) with a 
strength equal to unity, so that from the respective seed composition a subsequent r-process 
can occur. The third and seventh column contain the highest entropy (Sfi na i), for which 
seed nuclei can be formed. Beyond Sfi na i, no Fe-group seed nuclei are produced, and hence 
under such conditions an r-process is no longer possible. Obviously, F e -values larger than 
0.49 will exhibit a depletion of the third peak elements and the actinides Th and U since the 
maximum Y n /Y see( i ratios are smaller than 100. 

From Table H] one can also see that for a given expansion speed, Sf ina i changes and increases 
as Y e decreases. This is because the bottle-neck reaction of Eq. [11] additionally depends on 
Y e , and - provided neutrons are available - it can proceed at densities smaller than those 
needed by the triple-a reaction since an uncharged particle is involved. 



THE R-PROCESS 



Starting with the conditions resulting from the charged-particle freeze-out presented in 
the previous sections, we have performed our r-process calculations with the code discussed 
in Section 2. In addition to the original code version of iFreiburghaus et al.l (119991 ) . we also 
have introduced temperature-dependent neutron capture rates. We have investigated the 
impact of different mass models and /3-decay properties on the r-process duration and the 
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final abundances after /3-decay back to stability, and also have examined the role of /3-delayed 
neutrons. Special emphasis was put on a detailed understanding when a chemical equilibrium 
between neutron captures and photo-disintegrations, i.e., an (n, 7) — (7, n)-equilibrium is 
fulfilled, when the freeze-out from such conditions occurs and how neutron captures of the 
remaining free neutrons and /3-delayed neutrons can affect the final abundance pattern. 



4.1. General Results 

Here we present a first survey of detailed network calculations for a variety of entropies, 
followed until the total consumption of neutrons. According to the entries in Table HI en- 
tropies of 175, 195, 236, 270 and 280 k B /baryon result in neutron to seed ratios of about 
23, 34, 66, 107 and 121. The first two components lead to abundance patterns with maxima 
between A=100-130. The third component populates the rare earth elements (hereafter, 
REE) in the mass range around A=162. The fourth component forms the abundances of 
the A=195 peak, and the fifth component proceeds to even heavier isotopes, also populating 
the actinides. These results are displayed in Figs. [6] to [TO], making use of five different mass 
models/formulae. Fission is here implemented in a simplified way, assuming complete and 
symmetric fission for nuclei with A > 260. With a superposition of different entropy yields 
under the assumption that the ejected neutrino wind elements per equidistant entropy inter- 
val have equal volumes, final isotopic r-process abundances are obtained, which are displayed 
in Figs. [21] to [28j Deficiencies to the N r> Q will be discussed in more detail later in section 
5.2. 



4.2. Neutron Freeze-out 

While the r-process proceeds, the neutron abundance Y n decreases and the amount of 
the newly formed fresh r-process material increases as function of time. Since neutron cap- 
tures destroy the seed nuclei, we replace the seed abundance Y see d by the r-process material 
abundance Y r . With this, we define the neutron freeze-out as the instant when 

m < l - (13 > 

and we label the corresponding time as tf reeze . After tf reeze , the remaining neutrons and those 
produced by /^delayed neutron emission will not change the final abundances dramatically. 
Figs. 161 to [TU1 show for 4 mass models and one mass formula the freeze-out times tf reeze for 
five selected entropies, i.e., S = 175,195,236,270 and 280, and the final abundance curves 
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after /3-decay back to stability. Table |5] shows a comparison between two specifically selected 
mass models, i.e., FRDM and ETFSI-Q, for the same astrophysical conditions (V exp , Y e and 
S). Obviously, the nuclear physics input seems to have a significant impact on the r-process 
abundances. For instance, with the FRDM mass model a longer process duration (by about 
a factor 1.5) is needed to produce the second (A ~ 130) and third (A ~ 195) r-process peaks 
than for ETFSI-Q. Similar differences in the time scales are also observed for ETFSI-1, 
DUFLO-ZUKER and HFB-17. These differences are mainly due to the strong N = 82 shell 
closure of the FRDM mass model which in addition leads to an underproduction of the REE 
mass region. 

The weaker N = 126 shell closure of FRDM causes an underproduction of Pb and an earlier 
onset of fission cycling than the other mass models. Hereafter, when speaking about freeze- 
out without any further specifications, the neutron freeze-out (Y n /Y r < 1) is meant. 



4.3. Chemical Freeze-out 

In case of a chemical equilibrium between neutron captures and reverse photodisinte- 
grations, or a so-called (n, 7) — (7, n)-equilibrium between neighboring nuclei (Z, A) and 
(Z,A + 1) within an isotopic chain, the abundance ratio of these two nuclei is given by the 
nuclear Saha-equation 

Y(Z,A + 1) _ n n G(Z,A + l) / A+l \ 3/2 / 2nh 2 \ 3/2 f -S n (Z,A+l) \ 
Y(Z,A) 2 G(Z,A) \ A J \m u k B T ) 6XP V k B T ' ' 

The neutron separation energy S n , the nuclear partition function G and the mass number 
A represent the nuclear-physics input, whereas the temperature T and the neutron number 
density n n represent the conditions of the astrophysical environment. In the following, we 
want to analyze our r-process calculations, testing if in high-entropy expansions an (n, 7) — 
(7, n)- equilibrium is achieved at all, how long it prevails and at which time it breaks down. 
In addition, we will monitor the effect of abundance changes after the freeze-out from such 
an equilibrium via final captures of remaining free neutrons and those from /^-delayed netron 
emission. If we express the product k B T in units of 9 , Eq. [13] can be rewritten as 

nM±i) _ 8 . 410 -35 nn g^+l) ( ^±l) 3 ' 2 r -3/, exp (_ H-604 W A + l) 



Y(Z,A) G(Z,A) V A 

From Eq. [TS] we then obtain 



(15) 



S£^' n \Z,A+l) = -^ l \n 



12iqU Y(Z,A+1) G(Z,A) fA + iy 3/2 T* /2 



Y(Z,A) G(Z,A + 1) \ A J n. 



6) 
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S n (Z,A + 1) can also be calculated via 

S r n eal (Z, A + 1) = mf/ +1 - (mf/ + m n J, (17) 

where m^. +1 , m^ A and m^. are the mass excesses of the nuclei (Z, A + 1), (Z, A) and the 
neutron, respectively. We define the freeze-out with regard to the (n, 7) — (7, n)-equilibrium 
as the instant, when the ratio of Si™' 7 7 ' n ' ) and S^ eal of the nuclei (Z, A + 1) next to those 
lying in the r-process path (Z, A) (which represent the abundance maxima in each isotopic 
chain computed via relations given in Eqs. [TH and [T7] diverge from unity, and we label the 
corresponding time as t c h em . 

Figs. [TT1 to [TBI show the validity and duration of the (n, 7) — (7, ra)-equilibrium as a function of 
time for the entropies which synthesize the second, the REE pygmy and the third r-process 
abundance peaks. Obviously, the second peak is formed under a full chemical equilibrium. 
The ratios of the real neutron separation energies and those predicted by the Saha-equation 
are about unity for the whole Z-range (blue color). The break-out from the chemical equi- 
librium happens 60 ms later, i.e., clearly after the neutron freeze-out. Therefore, for the 
A ~ 130 r-process peak one has 

t freeze ^ tchem ■ 

The broad REE pygmy peak with its abundance maximum around A ~ 162 is formed at the 
limit of a chemical equilibrium, where 

t freeze ~ ^chem ■ 

However, the A ~ 195 peak is no longer built under chemical equilibrium conditions. This 
means 

t freeze ^* ^chem • 

Therefore, the morphology of the rapid neutron-capture process in the neutron-rich HEW 
can be depicted combination of three types: 

1. a "hot" r-process, which proceeds and ends already at relatively high temperatures 
(Tg rs 0.8), maintaining a full chemical equilibrium until the A ~ 130 peak is formed. 

2. a "hybrid" r-process, which proceeds at high temperatures and ends at low temper- 
atures (T g rj 0.5), maintaining a partial chemical equilibrium until the REE pygmy 
peak is formed. 

3. a "cold" r-process, which proceeds at high temperatures and ends at even lower tem- 
peratures (T^_R£_0 1 35jwi tlijw cherm until the A ~ 195 peak is formed 



see also (I Ar cones et al .1120071 ; IWanajo 1 120071 )). 



Thus, for the latter two r-process types the neutron capture rates will play a significant role 
for the mass region beyond A ~ 140. 
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4.4. Dynamical Freeze-out 



The dynamical freeze-out occurs when the neutron depletion becomes inefficient, and 
thus the timescale r n is larger than the hydrodynamical timescale Thyd, i-e., 



Y 



> 



P(t) 



P(t) 



18) 



tdyn — t freeze is the time interval 



We label the corresponding time with td yn - Then, At 
when the last free neutrons are captured and /3-delayed neutrons are recaptured during the 
decay back to stabillity. Depending on the capability of the radioactive progenitor nuclei to 
emit /3-delayed neutrons after tf reeze , the freeze-out times td yn and tf reeze may be identical 
or can differ from each other. Beyond td yn the neutron number density n n drops to values 
< 10 18 n/cm 3 , thus representing the signature for the end of the r-process. 



4.5. The /3-Flow Equilibrium 



During the r-process the abundance flow from each isotopic chain to the next is governed 
by /3-decays. We can define a total abundance in each isotopic chain Y(Z) = ^2 A Y(Z, A), 
and each Y(Z,A) can be expressed as Y(Z,A) = Y(Z)P(Z, A), where P(Z,A) represents 
the population coefficient of the nucleus (Z,A). In case of a steady flow of /3-decays between 
isotopic chains, in addition to an (n, 7) — (7, n)-equilibrium, we also have 



Y(Z, A)\ P (Z, A) = Y(Z) P(Z, A)\p(Z, A) 



Y(Z)Xf(Z) 



const. 



(19) 



In other words, each Z-chain can be treated as an effective nucleus with a /3-decay rate 
\ e ff. Under such conditions the assumption of an abundance Y(Z m i n ) at a minimum Z- 
value would be sufficent to predict the whole set of abundances as a function of A, provided 
that an (n, 7) — (7, n)-equilibrium and a steady flow has been reached. Already in the 
past, various groups have discussed the possibility of either a global r-process steady- flow 
equilibriu m or local (n, 7) — (7, n) and /3-flow equ i libria in between neutron shell closures 
' see, e .g. , ( jBurbidge et al.l 119571 ; ISeeger et al.l Il965l ; iHillebrandt et al.l Il976l ; ICameron et al. 



1983aJJb()), assuming different astrophysical sites and different theoretical nuclear-physics 
input. First experimental evidence for a combination of an ( n, 7) — (7, n) - equilibrium and a 
local steady flow at N=50 and N=82 came from the work of kratz et all Jl986i Il988l Il993h 
in their site-independen t approach. T his result that the /3-flow is not global was shortly 



afterwards confirmed by iMeyer I ( 119931 ) for the specific site of a high-entropy hot bubble. In 
the present paper, we extend the above steady-flow studies and show under which entropy 
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and temperature conditions which type of equilibrium is reached and when it breaks down. 
In Figs. [T4iri6l we show the quantity Y(Z)X e J (Z) (color-coded) as a function of time for 
entropies of S = 195, 236 and 270 which produce the second peak, the REE pygmy peak 
and the third r-process peak regions, respectively. Those entropies were selected because 
they indicate best how the /3-flow proceeds through both shell closures (N=82 and 126) 
and the REE region between them. As a function of time, the following conclusions can be 
drawn: 



Before approaching the freeze-out, a /3-flow equilibrium is established for isotopic chains 
with Z<43, i.e., where the abundance distribution does not populate the N=82 shell 
closure. 

After reaching the freeze-out, another late /3-flow proceeds through the most populated 
Z-chains: (a) Z=46-49 for the A~130 mass region, (b) Z=50-63 for the REE pygmy 
peak region and (c) Z=64-70 for the A~195 mass region. 



Thus, we confirm the earlier conclusions (see, e.g., (IKratz et al.l Il993l ; iMeyer I Il993 



Thielemann et al.l Il994l ; iPfeiffer et al.l l2001al )). that in the r-process a global /3-flow equi- 
librium does not occur, but local equilibria are established in between neutron shell closures 
where short /3-decay half-lives are encountred. These equilibria are essential for a succesful 
reproduction of the N riQ distribution (see Chapter 5). 



4.6. The Progenitors of 130 Te, 162 Dy and 195 Pt, and the Role of /3-Delayed 

Neutrons 

130 Te, 162 Dy and 195 Pt represent the maxima of the A~130, the REE and the A~195 
r-process abundance peaks. While it is generally accepted that the two sharp main peaks 
(both with HWFM widths of about 15 m.u.) are due to the retarded r-matter flow at 
the N=82 and N=126 shell closures, the origin of the broad structured REE pygmy peak 
distribution (with a total width of about 50 m.u.) as either being a signature of nuclear de- 
formation or as representing the result of a late abundance enhancement by fission cycling, 
still seems to be under debate. In the following, on the basis of Table [6] and Figs. [1711281 
we will discuss in some detail the formation of the three r-abundance peaks. Table [6] shows 
the isobaric progenitors of the three peak maxima at A=130, 162 and 195, respectively, for 
the three freeze-out types at different times, temperatures and neutron densities discussed 
above. Figs. [TTUT9l indicate the effects of the late recapture of /3-delayed neutrons on the 
final r-abundances. Figs. 12011231 show the calculated abundance heights and widths as a func- 
tion of entropy of the above three peak isotopes for the two selected mass models ETFSI-Q 
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and FRDM. And finally, Figs. I21H281 present overall fits to the N rj0 distribution for the five 
different mass models chosen in the present study. 

In the historical WP approximation of the r-process at the commonly assumed freeze-out 
around Tg~ 1, the main isobaric progenitor of stable 130 Te is the N=82 isotope 130 Cd, for 
which in the meantime th e most important nuclear- physics properties (T1/2, P n and Qg ) 
are experimentally known (IKratz et al.lll986t lHannawald et al.ll2000t iDillmann et al.ll2003l ). 
In a fully dynamical r-process model like our present HEW study, however, the situation is 
more complicated. In the following, we will discuss in detail the formation of 130 Cd during 
the freeze-out phases (see Table [6] and Figs. [11], [14] and [TT] For this purpose, we consider 
here the (representative) conditions of an entropy of S — 195, for which the highest partial 
130 Te abundance is obtained. 

Under these conditions, the inital r-matter flux at a neutron density of n n ~ 10 27 occurs 
further away from the /3-stability line than in the later freeze-out phases. As in the classical 
WP a pproach, the r-process in the A~130 region still starts to "climb up the N=82 stair- 
case" (IBurbidge et al.lll957l ). however only up to Z=45 127 Rh. Already at Z=46 128 Pd, the 
r-process begins to break out from the magic shell, with comparable time-scales for the two 
competing reactions of further neutron capture and /3-decay. In the Z=46 isotopic chain, the 
r-process does produce some N=82 128 Pd, but proceeds up to N=84 130 Pd (with the highest 
isotopic abundance) and N=86 132 Pd. Similarly, also for Z=47 the r-process partly breaks 
out from N=82, still producing 129 Ag, but again continuing up to N=84 131 Ag and N=86 
133 Ag with comparable initial abundances. In this early stage, for the classical Z=48 130 Cd 
waiting point, neutron capture out of N=82 still dominates over /5-decay (r„ 7 /r^ ~ 0.1) thus 
forming only a small fraction of its total abundance. Here, the highest initial isotopic yield 
is obtained for N=86 134 Cd. Finally, the classical N=82, Z=49 "break-out" isotope 131 In is 
passed quickly (r n7 /rg ~ 0.02), and N=86 135 In is formed with the highest initial isotopic 
yield. 

When focussing in the following on A=130, we recognize that several isotopes lying beyond 
N=82 with high initial abundances act as /3-decay and /3-delayed In to 3n progenitors of 
neutron magic 130 Cd. As mentioned above, the most important isobaric r-process nuclide 
now is N=84 130 Pd, which contains about 60% of the total abundances of the early progen- 
itors of 130 Cd. However, due to its high /3-delayed neutron branching ratios (P in ~ 73%; 
P 2n ~ 18%), and the total P xn ~ 24% of the daughter 130 Ag, only about 6% of the initial 
130 Pd abundance lead to 130 Cd. The next important initial /3-xn progenitors of 130 Cd are 
132 Pd and 133 Ag, both with about 20% progenitor abundance. With these different decay 
branches, during the early neutron and chemical freeze-out phases (see Table [6]), only a part 
of the total 130 Cd abundance is produced. Another sizeable fraction comes from neutron 
captures of the remaining "free" neutrons. Here again, the most important contribution 
comes from the initial 130 Pd abundance via its stong /3-delayed In decay branch to 129 Ag. 
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When considering the different neutron-capture cross sections for N=82 129 Ag and its /3- 
decay daughter N=81 129 Cd (about a factor 55 in favor of 129 Cd) under the relevant time, 
temperature and neutron-density conditions at the end of the chemical freeze-out (see Ta- 
ble E]), it is clear that neutron capture predominantly occurs after 129 Ag /3-decay to 129 Cd 
into N=82 130 Cd. The respective production modes of 130 Cd are reflected by the abundance 
ratios of Y( 130 Cd)/Y( 129 Ag) at the beginning and at the end of the chemical equilibrium. 
In this freeze-out phase, within about 240 ms, this ratio changes by about a factor 20 from 
0.14 to 2.7. These early abundances are further modulated during /3-decay back to stability 
by the capture of the remaining "free" neutrons and the /3-delayed neutron recapture. The 
resulting abundance shifts in the A~130 peak region from the late recapture of previously 
emitted /3-delayed neutrons are shown in Fig. [T71 We see that significant effects occur in the 
rising wing of the peak. At the top, the abundance of 129 Xe is reduced by a factor of 2.5, 
whereas the 130 Te abund ance increases by about a factor 2.2. 

Early debates (see, e.g., (IBurbidge et al.lll957l ; ICameron 1119571 )) were related to the possi- 
ble orogin of the REE pygmy peak, either due to mass-asymmetric fission cycling from the 
trans-actinide region or fr om the nuclear sha pes of the REE progenitor isotopes. In their 
site- independent approach iKratz et all ( 119931 ) had already recognized obvious deficiencies of 
the FRDM mass model in the shape-transition regions before and behind the N=82 and 
N=126 shell closures, which resulted i n an inadequat e REE pattern. More detailed network 
calculations were then performed by ISurman et al.l (119971 ). again using the FRDM mass 
model and a simplified set of "adapted" /3-decay halfe-lives. They concluded that the REE 
peak is due to a subtle interplay of nuclear deformation and /3-decay, not occuring in the 
steady-flow phase of the r-process. At about the same time, the conclusion of the above 
authors were in principle confirmed by r-abundance calculations of IKratz et al.l (119981 ). who 
compared the REE distri butions resulting from (i) the spherical mass mode l HFB/SkP of 



Dobaczewski et al.l ( 11984] ). and (ii) the deformed mass model ETFSI-Q of iPearson et al. 



(Il996f ). 

In the present paper, we use use five selected (deformed) mass models with different micro- 
scopic sophistication to reproduce the overall N riQ distribution. From Figs. I24"]|2"8l one can 
clearly see the different success of our dynamical network calcualtions in reproducing the 
overall shape of the REE pygmy peak. However, with the exception of FRDM , for the other 



four c ases we can consist e ntly c onfirm and strengthen the earlier conclusions of ISurman et al. 



( 119971 ) and IKratz et al.l ( 119981 ). that (at least for the HEW scenario) fission cycling seems 
to be unimportant. And, when considering the FRDM masses, it is highly questionable if 
its nuclear deficiencies around the magic neutron shells would be "repaired" more or less 
artificially by assuming strong fission cycling at extremely high entropies (e.g., S > 450 
k b /baryon). 

From Table [6] one can also see that the REE pygmy abundance peak is formed with the 
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"correct" shape and relative height in a continuous r- matter flow during feeze-out. It is in- 
teresting to note that in the early phase (begin of the chemical equilibrium at about 220 ms) 
the build-up of the whole REE region is still dominated by neutron captures, forming Z=54 
162 Xe and Z=56 162 Ba with the highest A=162 isotopic abundances. Within the following 
about 100 ms, when approaching the end of the chemical equilibrium, the matter flow has 
continued to the maximum A=162 abundance for Z=58 162 Ce. Now, up to this mass region 
/3-decay clearly dominates over further neutron capture. Only in the late freeze-out phase of 
the dynamical freeze-out (after about 480 ms), the continued matter flow has formed Z=59 
162 Pr and Z=60 162 Nd with their maximum yields. Under these low-temperature conditions 
(Tg ~ 0.3), the density of "free" neutrons has already dropped to n n ~ 10 17 so that by now 
the whole REE region between Z=54 and Z=62 is dominated by /3-decays. Finally, as is 
indicated in Fig. (TSJ during even later times of the decay back to stability, an additional 
gradual abundance shift to A=162 occurs, which originates from the recapture of previously 
emitted /3-delayed neutrons. Thus, the final shape of the REE pygmy peak is established 
only at very late non-equilibrium. 

The formation of the A~195 N rj0 peak, which is related to the N=126 shell closure, is similar 
to that of the A~130 (N=82) peak but not exactly the same. There are several differences 
with respect to both astrophysical as well as nuclear-physics parameters. Starting with a sim- 
ilar initial neutron density of n n ~ 10 27 , the higher entropy of S = 270 (compared to 5* = 195 
for the A~130 peak) leads to a higher neutron to seed abundance ratio of Y n /Y see d ~ 105 
(compared to Y n /Y see( i = 36 for 5* = 195). This results in an inital r-process path further 
away from /3-stability than in the A~130 region, involving extremely neutron-rich progeni- 
tor isotopes with on the average higher /3-delayed neutron branching ratios. In consequence, 
the description of the formation of the A~195 abundance peak and its modulation during 
the decay back to /3-stability is more complicated. Apart from the "bottle-neck" behavior 
of the N=126 shell closure, we have to face a detailed interplay between captures of "free" 
neutrons, /3-decays, emission of /3-delayed neutrons and their recapture. For some snapshots 
of the freeze-out during the first 750 ms, see Table [6] and Figs. [131 tHH an d [El And at the 
end, we will argue why the top of the A~195 N rj0 peak occurs with 195 Pt at an odd mass 
number, and not - as one would expect - at an even mass number. 

In principal similarity to the formation of the early A~130 N r peak, the initial r- matter 
flux enters the A~195 peak region at the "lighter" N=126 waiting points, producing small 
fractions of Z=62 188 Sm to Z=66 192 Dy. Already from Z=67 on upwards, the early r-process 
breaks through the closed shell. In this phase, in the whole mass region neutron capture 
strongly dominates over /3-decay. 

At the end of the chemical equilibrium freeze-out phase after about 420 ms (see Table E]), the 
abundances of N=126 191 Tb to 196 Yb have increased between one up to 7 orders of magni- 
tude, all approaching their maximum values during the first 750 ms of the freeze-out. For all 
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"lighter" N=126 isotopes up to Z=68 194 Er, now /9-decay dominates over neutron capture. 
Under these conditions, 194 Er has the highest N=126 isotopic abundance. For Z=69 195 Tm 
and Z=70 196 Yb neutron capture and /3-decay time scales become comparable; hence, the 
break-out from N=126 now has shifted to higher Z values. At the time when the dynamical 
freeze-out occurs (after about 750 ms; see Table [6]), the abundances of the lighter N=126 
isotopes up to Z=67 (Ho) have already decreased by more than an order of magnitude due 
to their /3-decays, whereas under these late conditions finally also the abundances of the 
"heaviest" N=126 isotopes from Z=71 (Lu) to Z=74 (W) have reached their maximum val- 
ues. Interestingly, the abundance of Z=68 194 Er has only decreased by about a factor 4, and 
the abundance of 195 Tm has remained constant. In this late dynamical freeze-out phase, for 
all N=126 nuclei /3-decay dominates over neutron capture. 

As mentioned above, at the end of this section we want to come back to the question why 
the top of the A~195 N r peak occurs at an odd mass number and not at the - as expected 
- even mass A=194 of the most abundant classical N=126 waiting-point isotope 194 Er. At 
the end of the chemical equilibrium and in the subsequent dynamical freeze-out phase, a 
subtle interplay of several /3-decay, /31n- to /33n-decay and neutron-capture channels occur, 
which lead to kind of reaction cycles. Recapture of /3-delayed neutrons does not yet play a 
significant role. Under these conditions, the maximum abundance in the peak is still 194 Er 
closely followed by 195 Tm. When neglecting /3-delayed neutron recapture, obviously this 
signature is held during the whole decay back to stability, resulting in a final abundance 
ratio of Y( 194 Pt)/Y( 195 Pt)~1.3. However, when allowing the recapture of previously emit- 
ted /3-delayed neutrons, which occurs relatively late during the /3-backdecay, a significant 
upwards abundance shift with Z occurs in the whole rising wing up to the top of the peak. 
Beyond A=195, for a few masses the reverse effect is observed. This is shown in Fig. [191 In 
our calculations, the final abundance ratio now changes to Y( 194 Pt)/Y( 195 Pt)~0.5. Hence, 
we again conclude that the recapture of /3-delayed neutrons has to be included in detailed 
r-process calculations. 



5. REPRODUCTION OF THE SOLAR-SYSTEM R-PROCESS 

ABUNDANCES 

5.1. The Entropy Weights 

In our attempt to not only understand the behavior of individual entropy components 
but rather to find an explanation of the solar system isotopic r-abundance residuals N r)0 ~ 
N — A s , we need to consider a recipe for the weights in such entropy superposit ions, similar 



to superposition of neutron number densities in classical r-process calculations (IKratz et al. 
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1993; 


Cowan et al. 


1999; 


Pfeiffer et al. 


2001b: 


Kratz et al. 


2007) 



two simple approaches, assuming the ejection of neutrino-wind elements per equidistant 
entropy interval with equal masses (dM(S)/dS = const) or with equal volume (dV(S)/dS = 
const), where the entropy spans from S = 5 to Sf ina i. This is a pure assumption not resulting 
from from dynamic explosion calculations; but we will show that this simple weighting leads 
to excellent fits to the solar r-process abundances between the rising wing of the A~130 peak 
and the Pb-Bi peak. 



5.1.1. Ejecta per equidistant entropy with equal volumes 

If the ejected neutrino- wind elements are equal-sized (hereafter, assumption 1), one 
obtains 

M(S) 

Since S ~ 1/p, one obtains for an entropy S 



const. 



M{S) S = M(S ref ) S ref = const, (20) 

where S re f represents an arbitrary reference entropy which we set equal to 5 kg/baryon, the 
smallest entropy in our network calulations. From Eq. [201 one infers that the entropy weight 

Hence, by integrating over different entropies, the accumulated abundance Y sum (Z,A) of a 
given nucleus [Z, A) is 

Y sum (Z,A)= V ^L Ys (Z,A). (22) 

5=5 ^ 

Therefore, the ejected mass fraction of heavy material per entropy is given by 

X heavy (S) = w(S)(l - X a ). 

Figs. M and ED show the calulated abundances of 130 Te, 162 Dy and 195 Pt function of 
entropy. The agreement of the peak heights of our simulations with the corresponding solar- 
system r-abundances is excellent for the mass models ETFSI-Q and FRDM. However, the 
latter mass model underproduces the REE mass region because of its very strong N=82-shell 
closure. 
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5.1.2. Ejecta per equidistant entropy with equal masses 

Here we show (as a different example) superpositions with equal amounts of ejected 
mass per entropy interval (hereafter, assumption 2). When the ejected masses per entropy 
are equal, the entropy weights will be the same, i.e., w(S) = 1. The accumulated abundance 
of a given nucleus then is 

S final 

Y sum (Z,A)= Ys(Z,A). (23) 

5=5 

Figs. [22] and [23] show the results for ETFSI-Q and FRDM. Both mass models show similar 
(but slightly worse) fits to the solar r-abundances and thus indicate that within the expected 
uncertainties both assumptions 1 and 2 yield good agreement. This might indicate what one 
would expect from realistic scenarios. For the attempt to reproduce the solar r-abundances, 
in the following we use the assumption 1. 



5.2. Fitting the Solar r-Process Abundances 



In the past, many attempts were undertaken to reproduce the solar r-process abun- 
dances. In the frame of the WP approximation, based on an (n, 7) — (7, n)-equilibrium 
with static temperatures and neutron number densities, it has been shown that only a su- 
perposition of a several of neutron number density components with weights given by the 
realistic heights of th e three N ri0 -peaks (Y( 8 Se)~ 24.3, Y( 130 T e)~ 1.69, Y( 195 Pt) ~ 0.457 



( iKappeler et al.lll989l )) can achieve this goal ( iKratz et al.lll9931 ). iFreiburghaus et al.l (119991 ) 
used an entropy-based superposition with an analytical fit function g(S) = xie~ X2S , where 
X\ and X2 represent two fit parameters. We follow the method discribed in this work by using 
our above assumption 1. Figs. [2H to US\ show the results for a selected constant expansion 
velocity of V exp = 7500 km/s and a specific electron abundance Y e = 0.45. It is surprising 
that this simple parametrisation of the HEW components provides a good fit to the overall 
solar r-abundances above A=110, especially with the mass models ETFSI-1 and ETFSI-Q. 
The mass formula of DUFLO-ZUKER slightly underproduces the A~195 and the Pb mass 
regions, and the new mass model HFB-17 shows deficiencies in the A~145 and A~175 mass 
regions where nuclear phase transitions occur. In these regions, the mass model prescriptions 
may still have systematic deficiencies. HFB-17 also underproduces the A~195 and the Pb 
mass regions. The mass model FRDM has defeciencies in producing the REE, the A~200 
and the Pb mass regions. 

The over-production of the mass region below A< 110 by the HEW is well-known and was 
discu ssed by several authors in the past (see, e.g., ( jWoosley et al.lll994l ; IFreiburghaus et al. 
19991 )). It has been referred to as a model deficiency due to the a-rich freeze-out component 
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where the r-process cannot start because the neutron to seed ratio is still smaller than unity. 
Fig. [291 shows the fit to the solar elemental r-abundances by using the mass model ETFSI-Q. 
It is quite evident from the this figure that the abundance over-production mentioned above 
concerns the elements between Sr and Ag. However, in our new understanding of the HEW 
nucleosynthesis one can avoid this disagreement by considering an additional superposition 
in terms of Y e . 

Since S increases and Y e decreases with time, it is reasonable to consider the ejecta of a 
core-collapse supernova explosion as a mixture of different S and Y e -components. Fig. I3"U1 
shows the mass fraction of the expected heavy ejecta (A> 4) for a selection of Y e - values as a 
function of S. Thereby, one can see that the weight of the ejecta increases as Y e decreases. 
Fig EH shows that the whole mass region from Sr (Z=38) up to U (Z=92) can be fitted by 
using 4 different Y e - values (0.498, 0,496, 0,490 and 0.482). Thereby, we have normalized 
the solar r-abundances so that Y(Nb) rj0 = F(Nb)y e=0 .498- We have chosen the element Nb 
because it can serve as normalization for the elemental and isotopic solar r-abundances since 
it has only one stable isotope. From Fig. [31] one can infer: 

• Sr, Y, Zr and Nb (Z=38-41) are formed under condition of Y e = 0.498 

• Mo and Ru (Z=42 and 44) are formed under conditions of Y e = 0.496 

• Rh, Pd and Ag (Z=45-47) are formed under further decreasing Y" e -conditions of Y e = 
0.490 

• Cd and the elements beyond (Z> 48) are formed under more neutron-rich conditions 
of Y e = 0.482 

Thus, the well-known over-production of the mass region from A= 88 to 110 in the HEW 
in Figs. [21] to [28] is due to the f act that only one Y e is used. A superposition of several 



Y e - values (as e.g., performed by iMeyer et al.l (119921 ) for a selected entropy), however in 



addition to the superposition of S- values used in our study, could help to resolve this problem. 



5.3. The Amount of the r-Process Ejecta 

Based on our assumptions concerning the entropy weights, which provide a good fit to 
the solar r-abundances beyond the A=110 mass region, we will estimate the amount of the 
(neutron-capture) r-process material which is ejected, 

M r « \ n Rl(t = 0) f Unal dS(l - AY a )p(S, t = 0). (24) 

,y ^initial 
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The range from SinUM to Sfi na i indicates the entropy interval which permits the production 
of r-process nuclei. Furthermore, our calculations show that the entropy-dependent r-process 
mass fraction X r = 1 — 4Y a can be well approximated by a function f(S) = XiS x ' 2 . Eq. |2H 
then reads 

M r « i 7T Rl (t = 0)xi J dS S X2 p(S, t = 0). (25) 
Replacing the term of the matter density in Eq. [25] by its term in Eq. El one obtains 

M r « 1.1 1O~ 3 M Xl J dS (26) 

in units of solar masses M & . The paramters xi and X2 vary with V^, and Y e . If we complete 
the integration in Eq. [261 one finally gets 

M r « 1.1 1(T 3 M (S^ - 5- tw ) . (27) 

Tables shows the results according to Eq. [27] for two expansion velocities. The amount 
of the estimated mass of the ejecta then depends only on the electron abundance Y e . For 
instance, at a given expansion speed & Y e = 0.45 will produce roughly 10 times more r-process 
material than 0.49. This is due to the higher entropy needed to start an r-process with the 
latter Y e - value. 



6. SUMMARY AND CONCLUSIONS 

We have performed large-scale parameterized dynamical network calculations in the 
context of an adiabatically expanding high-entropy wind as expected in core-collapse super- 
novae. We have used four mass models and a mass formula to test in great detail the validity 
of the (n, 7) — (7, n)- and the /9-flow equilibria and the time intervals during which they are 
fulfilled during an r-process. We have defined different kinds of freeze-out, (a) chemical, 
(b) neutron and (c) dynamical, that are related to (a) the break-out of the (n, 7) — (7, n)- 
equilibrium, (b) to the consumption of neutrons, i.e., Y n /Y heavy < 1, and (c) to the time at 
which the neutron depletion timescale becomes larger than the hydrodynamical timescale, 
i.e., Y n /Y n > \p(t)/p(t)\, respectively. In our attempt to provide the best possible fit to the 
solar r-process abundances, we have explored the consequences of two simple assumptions 
concerning the matter ejection from the newly born proto-neutron star in the HEW. The 
first assumption is based on ejection of neutrino-wind elements with equal sizes, and the 
second assumption is based on the ejection of neutrino-wind elements with equal masses per 
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equidistant entropy interval. Since the first assumption yielded exactly the same ratios of 
the solar r-abundance peaks ( 13O Te r . i0 / 162 Dy r . i Q w 16 and 13O Te n0 / 195 Pt ri ~ 3.5) we made 
use of it to determine the entropy weights, in order to superpose the entropy yields and 
to estimate the amount of the r-process material which can be ejected by the HEW of a 
core-collapse supernova explosion. 
Furthermore, we can draw the following conclusions: 



The entropies generated by the HEW exhibit maximum values after which the produc- 
tion of heavy seed nuclei via the bottle-neck nuclear reactions (Eqs. [10] and [TT]) fails. 
We have labeled those entropies with Sfi na i- The seed distributions beyond Sfi na i, con- 
sisting mainly of a-particles, neutrons, protons and traces of 12 C and 16 O, do not allow 
an r-process to proceed because of the bottle-neck regions at A = 5 and 8. Further- 
more, for a given expansion velocity of the matter in the HEW, the maximum entropies 
Sfi na i increase with decreasing Y e because the three-body bottle-neck nuclear reaction 
(Eq. [TT]) depends on the availability of neutrons in the wind. It can therefore proceed 
at lower densities than the triple-a-reaction since an uncharged particle is involved. 
An interesting observation is that (independent of Sfi na {) F e -values larger than 0.49 do 
not yield neutron to seed ratios large enough to synthesize the third r-process peak 
elements and Th and U. 

The mass region between Fe and Zr (depending somewhat on Y e ) which is historically 
thought to be the beginning of the "weak" r-process, is obviously formed by rapid 
charged-particle reactions in the HEW at low entropies. From about Nb on, our model 
predicts a smooth transition into a true rapid neutron-capture r-process. This seems, 
indee d, to be confirmed by recent astro-physic al and cosmochemical observations (see, 
e.g., fjKratz et al.l 120081 ; iFarouqi et al.l l2009al lbl) ) . Therefore, in this mass region the 
classical definition of the so-called Solar System (SS) r-process "residuals" as N — 
N 



s,0 



N r should no longer be applied. 



While in our calculations the lightest trans-Fe elements up to about Kr are underpro- 
duced by one to two orders of magnitude relative to the SS abundances, the region 
between Sr and Ag appears to be overproduced by about a factor 4 when using the 
full entropy range. These two kinds of differences relative to N r Q may be explained in 
the following way: (i) for the region between Fe and about Kr by m issing abundance 
contributions from the z/-p process (see, e.g., dFrohlich et al.l l2006ailbh_\ and/ or by ad- 



ditional contributions from a new type of rs-process (jPignatari et al.ll2008l ): and (ii) 
for the region between Sr and Ag by the assumption of a unique, constant electron 
abundance in the HEW when integrating over the entropy yields. 
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We also have shown by our HEW r-process studies that for electron abundances slightly 
below Y e = 0.50, e.g., for Y e = 0.498, only the mass region below the A~130 peak can 
be formed. The region of the classical "main" r-process up to the full 3rd peak requires 
somewhat more neutron-rich winds with Y e values in the range 0.48; and finally, a full r- 
process including sufficient abundances of the actinides Th and U can only be produced 
with even lower Y e values in the range below about 0.478. This is due to the fact that 
the amount of heavy ejecta as a function of S decreases with increasing Y e . 

The SS r-process residuals beyond A~110 can be well fitted by assuming that the HEW 
of core-collapse supernova explosions will eject a series of neutrino wind elements per 
equidistant entropy interval. Thereby, the relative heights of the A~130 and A~195 
major peaks as well as the intermediate REE pygmy peak in the r-abundance fits repro- 
duce the N ri0 distributions quite well for a variety of microscopic global mass models. 
Depending on Y e , the total r-process masses lie between bet ween 10~ 6 and 10~ 4 Mm, in 



good agreement with Galactic chem ical evolution studies ( ITruran &: Cameron 111971 



Hillebrandt 


1978; 


Cowan et al. 


1991) 



eters have sizeable effects on the process duration to overcome the major bottle neck 
in the r-process matter flow at the N=82 shell thereby forming the A~130 abundance 
peak, as well as the total r-process duration to produce sufficient amounts of Th and 
U. For example, by using the same astrophysical parameters (Y e = 0.45, S < 280 and 
Vexp — 7500 km/s) for the quenched mass model ETFSI-Q the full r-process up to the 
actinides is about a factor three "faster" than for the older, unquenched mass model 
FRDM. Therefore, for any type of modern r-process calculation in whatever astropy- 
hsical scenario, a sophisticated nuclear-physics input, in particular around the magic 
neutron shells, remains crucial. 

The morphology of the r-process in the HEW can be depicted as: 

1. A "hot" r-process (in full chemical equilibrium), which produces the mass region 
up to the A~130 r-process peak; 

2. a "hybrid" r-process (in partial chemical equilibrium), which is responsible for 
producing the REE region; and 

3. a "cold" r-process (out of chemical equilibrium), which forms the A~195 peak 
and the actinides. 

The REE pygmy peak of the SS r-process abundances in between the two major 
peaks at A~130 (caused by the N=82 shell closure) and at A~195 (caused by the 
N=126 magic shell) is formed in rather small entropy ranges for the more recent micro- 
scopic mass models, with the exception of the still frequently used older macroscopic- 



-27- 



microscopic model FRDM. With this, our results clearly show that with an appropriate 
choice of the nuclear-physics input, fission cycling is not neceessary to reproduce the 
correct shape of the REE pattern. 

• Without considering effects from the emission and recapture of /3-delayed neutrons, 
the two main r- abundance peaks with their tops at A=130 and A=195 would occur at 
A=129 and A=194. The late recapture of /3-delayed neutrons is the main reason for the 
gradual shifts from A=129 to the classical major N=82 "waiting-point" isotope 130 Cd 
in a quasi-chemical equilibrium, and the respective shift from A=194 to the classical 
N=126 "waiting-point" nuclide 195 Tm under non-equilibrium conditions. 

• The /3-flow equilibrium is confirmed not to occur globally during the whole r-process. 
Only local equilibria are achieved, partly at the magic neutron shells and in between 
the different shell closures. 
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Fig. 1. — Time evolution of temperature T and density p in the expanding hot bubble. The 
initial values of T and p are normalized to 1. The freeze-out of the charged-particle reactions 
occurs when the initial temperature in units of 10 9 drops from Tg = 9 to 3, and the freeze-out 
of a subsequent r-process occurs when Y n /Y r ^^>4) < 1. 




Fig. 2. — Production and depletion of the a-particles as a function of decreasing temperature. 
The maximum mass fraction of the a-particles which can be reached is X max = 2 Y e . The 
depletion of the a-particles depends both on the entropy S and the electron abundance Y e . 
Beyond a maximum entropy Sfi na i, the production of any heavy seed is not possible because 
of the very small density. Hence, under such conditions an r-process is no longer possible. 



-33 - 



x 

w 

<~ 

X" 
ll_ 

> 
> 

<d" 

Q 

c 
cd 

"O 

c 

3 
< 



High-Entropy Wind: V exp = 7500 km/s, Y e = 0.45 



0.1 



0.01 



Neutrons 
ct-Particles 



Seed 



0.001 r 
0.0001 — 



50 100 150 200 250 300 350 
Entropy, S(k B /baryon) 



Fig. 3. — Remaining abundances of a-particles, heavy "seed" nuclei and free neutrons after 
the charged-particle freeze-out as a function of the entropy for a selected electron abundance 
and expansion velocity. The r-process starts when the neutron to seed ratio becomes larger 
than unity, here at S « 100. 



High-Entropy Wind: Y e = 0.45, Y n /Y seed = 50 
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Fig. 4. — The seed nuclei distributions for an arbitrary constant neutron to seed ratio of 
Y n /Y see d = 50, an electron abundance Y e = 0.45 and for 4 different expansion velocities of 
the hot bubble, ranging from V exp = 3750 up to 30000 km/s. The corresponding entropies, 
ranging from S = 140 up to 270 kg/baryon are chosen such that the ratio mentioned above 
remains constant. The resulting seed distributions look very similar and are robust. The 
corresponding mass fractions of the a-particles X a are equal. 
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Fig. 5. — The seed nuclei distributions for an arbitrary constant neutron to seed ratio of 
Y n /Y see d = 50, an expansion speed V exp = 7500 km/s and 4 different electron abundances 
ranging from Y e = 0.45 up to 0.49. The corresponding entropies are chosen such that the 
ratio mentioned above remains constant. The resulting seed distributions are not robust. The 
corresponding mass fraction of the a-particles decreases with decreasing electron abundance. 
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Fig. 6. — A choice of five entropies which sy nthesize the mass region beyond A=110 by using 
the mass model FRDM (jMoller et al.lll995l ) and the corresponding /3-decay properties and 
neutron-capture rates. The solar r-abundances (black dots with error bars) are scaled such 
that the value of 13O Te rj0 coincides with the maximum of the S = 195 curve at A=130. The 
process durations indicate the time between the freeze-out of the charged-particle reactions 
and the r-process. Note the depletion ot the REE mass region. 
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HEW, ETFSI-1 , V exp = 7500 km/s, Y e = 0.45 
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Fig. 7. — A choice of five entropies whi ch synthesize the mas s region beyond A=110 by using 
the unquenched mass model ETFSI-1 (lAboussir et al.lll995l ) and the corresponding /3-decay 
properties and neutron-capture rates. See caption of Fig. [6] and text for further details. 



x 

w 

<f 

X" 
ll_ 

> 

CD 
o 
c 
CO 



< 



0.1 

0.01 
0.001 
0.0001 
1e-05 
1e-06 
1e-07 



HEW, ETFSI-Q, V exp = 7500 km/s, Y e = 0.45 

7 — ■ pz nlr0 py Process Duration- 

is 1 75, 086 ms 

f * ■, 195, 138 ms 

236, 253 ms 

270,398 ms 
280, 468 ms 




60 80 100 120 140 160 180 200 220 240 
Mass number, A 



Fig. 8. — A choice of five entropies whi ch synthesize the ma ss region beyond A=110 by using 
the quenched mass model ETFSI-Q (IPearson et al.l Il9961 ) and the corresponding /3-decay 
properties and neutron-capture rates. See caption of Fig. [6] and text for further details. 



-36 - 



HEW, Duflo-Zuker, V exp = 7500 km/s, Y e = 0.45 
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Fig. 9. — A choice of five entropies whi ch synthesize the mas s region beyond A=110 by using 
the mass formula of DUFLO-ZUKER flDuflo fc Zuker Ill996l ) and the /3-decay properties and 
neutron-capture rates calculated based on ETFSI-Q. See caption of Fig. |6]and text for further 
details. 
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Fig. 10. — A choice of five entr opies which synthesi ze the mass region beyond A=110 by using 
the new mass model HFB-17 (IGoriely et al.l 120091 ) and the /3-decay properties and neutron- 
capture rates calculated based on ETFSI-Q. See caption of Fig. [6]and text for further details. 
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The (n,Y)-(Y,n)-Equilibrium: S n real /S n (n ' YHy ' n) 
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Fig. 11. — The validity and endurance of the (n, 7) — (7, n)-equilibrium for the A=130 peak 
formed by S=195 kg/baryon. Color-coded is the ratio of the neutron separation energies 
S^ eal (Z, A+l) of the right neighbor of the nucleus (Z, A) with the maximum abundance in the 
corresponding Z-chain, calulated by using the ETFSI-Q mass model and ^ n ' 7)_(7,n) (Z, A+l) 
predicted by the nuclear Saha-equation in case of an existing chemical equilibrium between 
the nuclei (Z,A) and (Z,A + 1). The equilibrium holds when S^ eal /Sn ~ 1- The 

vertical solid line shows the time of the r-process freeze-out. See text for further details. 
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Fig. 12. — The validity and the endurance of the (n, 7) — (7, n)-equilibrium for the A=162 
REE pygmy peak formed by S = 236 k B /baryon. See caption of Fig. [TTJfor further details. 
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Fig. 13. — The validity and the endurance of the (n, 7) — (7, n)-equilibrium for the A=195 
peak formed by S = 270 kg/baryon. See caption of Fig. [LTJfor further details. 



The p-Flow Equilibrium: A, eff (Z) Y(Z) 
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Fig. 14. — The validity and endurance of the /3-flow equilibrium for the A=130 peak formed 
by S = 195 kg/baryon as a function of time. Color-coded is the quantity X e Q(Z)Y(Z) of 
each Z-chain. The vertical line shows the time of the r-process freeze-out. Here, a global 
/3-flow equilibrium does not occur, but local equilibria do appear between the populated 
odd-Z and even-Z chains before the r-process freezes out. After the freeze-out the /3-flow 
equilibrium proceeds steadily to higher Z-chains (Z=47 and 48) around t ~ 150 ms, and 
proceeds further to Z=50 before it completely breaks out. 
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The p-Flow Equilibrium: A. eff (Z) Y(Z) 
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Fig. 15. — The validity and endurance of the /3-flow equilibrium for the REE pygmy peak 
formed by S = 236 ke/baryon as a function of time. Prior to the freeze-out, a wide-range 
equilibrium occurs between the Sn (Z=50) and Nd (Z=60) Z-chains when the matter flow 
passes the In (Z=49) Z-chain at t ~ 200 ms. After the freeze-out, the equilibrium lasts about 
50 ms between the most populated Z-chains in the REE region. See caption of Fig. HH and 
text for more details. 
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Fig. 16. — The validity and endurance of The /3-flow equilibrium for the A=195 peak formed 
by S = 270 kg/baryon as a function of time. Prior to the freeze-out, a strong local equi- 
librium occurs between the Sm (Z=62) and Er (Z=68) Z-chains when the mater flow passes 
the Pm (Z=61) Z-chain at t « 250 ms. After the freeze-out, a weak equilibrium lasts for 
about 300 ms between the Er (Z=68), Tm (Z=69) and Yb (Z=70) Z-chains. See caption of 
Fig. dH and text for more details. 
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ETFSI-Q, V exp = 7500 km/s, Y e = 0.45, S= 195 
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Fig. 17. — The upper panel shows two abundance distributions for the same astropysical 
parameters (S, Y e , V exp ) with and without (—j3) recapturing of the previously emitted 
/3dns. Their recapture causes the maximum abundance to shift from A=129 to A=130. The 
lower panel shows the ratio of the isotopic abundances for both cases, without and with 
/3dns. See text for further details. 

ETFSI-Q, V = 7500 km/s, Y e = 0.45, S= 236 
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Fig. 18. — The upper panel shows two abundance distributions for the same astropysical 
parameters (S, Y e , V exp ) with and without (— j5) recapturing of the previously emitted 
/3dns. Their recapture causes the gradual shift of the maximum abundance to A~162. The 
lower panel shows the ratio of the isotopic abundances for both cases, without and with 
/3dns. See text for further details. 
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Fig. 19. — The upper panel shows two abundance distributions for the same astropysical 
parameters (S, Y e , V exp ) with and without {—0) recapturing of the previously emitted 
/3dns. Their recapture causes the gradual shift of the maximum abundance from A=194 to 
A=195. The lower panel shows the ratio of the isotopic abundances for both cases, without 
and with /3dns. See text for further details. 
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Fig. 20. — The calculated abundances of 130 Te, 162 Dy and 195 Pt as a function of the entropy 
S by using the quenched mass model ETFSI-Q. The three curves are scaled with the same 
factor to normalize the maximum of the 130 Te curve to the abundance of 13O Te rj 0. Note the 
good agreement of the the relative heights of the peak maxima with the corresponding solar 
r-abundances in case of a continuous ejection of neutrino wind elements with equal volumes 
per equidistant entropy interval (V(S) = const). 
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FRDM, V = 7500 km/s, Y e = 0.45 
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Fig. 21. — See Fig. [20j The mass model used here is FRDM. Again, note the good agreement 
of the maximum of the 195 Pt curve with its corresponding solar r-abundance. The underpro- 
duction of the REE mass region is due to the strong N=82-shell closure of the mass model 
used. 
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Fig. 22. — See Fig. [201 but with the ejection of neutrino-wind elements with equal masses 
per equidistant entropy interval (M(S) = const). Note the slight disagreement of the top of 
the 162 Dy and 195 Pt curves with their corresponding solar r-abundances. 
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FRDM, V = 7500 km/s, Y e = 0.45 
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Fig. 23. — See Fig. [20j but with the ejection of neutrino-wind elements with equal masses 
per equidistant entropy interval (M(S) = const) and by using the FRDM masses. Note the 
strong underproduction of 162 Dy and the slight disagreement of the top of the 195 Pt curve 
with its corresponding solar r-abundance. 
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Fig. 24. — Superposition of entropies from S = 5 up to the maximum entropy Sf ina i = 355 
using the mass model FRDM. The overall fit of the mass region from A=110 to 238 is 
good, except for the REE mass region which is strongly underproduced by a factor of 3. 
Furthermore, the left wing of the A=195 peak is displaced and an abundance trough occurs 
at A^205 which leads to a strong underproduction of the Pb mass region. 
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Fig. 25. — Superposition of entropies from S = 5 up to the maximum entropy Sf ina i = 355 
using the unquenchend mass model ETFSI-1. The overall fit of the mass region from A=110 
to 238 is excellent. 
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Fig. 26. — Superposition of entropies from S = 5 up to the maximum entropy Sf ina i = 355 
using the quenchend mass model ETFSI-Q. The overall fit of the mass region from A=110 to 
238 is good. However, the mass region between A=140 and 160 is somewhat overproduced. 
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Fig. 27. — Superposition of entropies from S = 5 up to the maximum entropy Sf ina i = 
355 using the mass formula of DUFLO-ZUKER. The overall fit provides the best overall 
agreement for the REE mass region, but the third r-process peak and the Pb mass region 
are underproduced. The abundance trough at A=205 is less pronounced than in the case of 
FRDM. 




Fig. 28. — Superposition of entropies from S = 5 up to the maximum entropy Sf ina i = 355 
using the mass model HFB-17. Note the slight underabundances in the nuclear shape- 
transition regions at A~145 and A~175. The third r-process peak and the Pb mass region 
are underproduced. The abundance trough at A=205 is more pronounced than in the case 
of DUFLO-ZUKER. 
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a- and r-Process Yields, Y e = 0.450, V exp = 7500 
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Fig. 29. — Solar-system elemental abundances obtained with the same method discribed in 
Fig. |25j Significant deviations of the fit from the solar r-abundances occur at the elements 
below Cd (Z<48). For further discussion, see text. 




Fig. 30. — Mass fractions of the expected heavy ejecta (A>4) as a function of entropy 5* and 
electron abundance Y e in case of ejection of neutrino-wind elements per equidistant entropy 
interval. The mass fraction of the heavy ejecta decreases with increasing S and Y e . 
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a- and r-Process Yields, ETFSI-Q, V exD = 7500 
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Fig. 31. — Solar-system elemental yields are compared to the results of 4 electron abundances 
ranging from 0.498 down to 0.482. For each Y e the superposition of the entropies spans from 
S = 5 to the maximum entropy Sf ina i(Y e ). The elemental solar r-abundances are scaled such 
that y r Q (Nb) = 3^0.498 (Nb). The element Nb is chosen as reference because it has only one 
stable isotope and therefore can also be used as reference for the isotopic solar r-abundances. 
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Table 1: Dimension of the full reaction network used until charged-particle freeze-out. The 
particle unstable isotopes 5 He, 8 Be and 9 B are included implicitely as intermediate nuclei 
in (formal) two-body representations of three-body reactions. The network contains 1989 
nuclei and the neutron. 



Element 


z 


A • 
^mm 


f\ 1 VI Q V 


Element 


z 


A • 
^mm 


1 MlloLX 


H 


1 


1 


3 


Cr 


24 


38 


86 


He 


2 


3 


6 


Mn 


25 


40 


89 


Li 


3 


6 


9 


Fe 


26 


42 


92 


Be 


4 


7 


12 


Co 


27 


44 


96 


B 


5 


8 


14 


Ni 


28 


46 


99 


C 


6 


9 


18 


Cu 


29 


48 


102 


N 


7 


11 


21 


Zn 


30 


51 


105 





8 


13 


22 


Ga 


31 


53 


108 


F 


9 


14 


26 


Ge 


32 


55 


112 


Ne 


10 


15 


41 


As 


33 


57 


115 


Na 


11 


17 


44 


Se 


34 


59 


118 


Mg 


12 


19 


47 


Br 


35 


61 


121 


Al 


13 


21 


51 


Kr 


36 


63 


124 


Si 


14 


22 


54 


Rb 


37 


66 


128 


P 


15 


23 


57 


Sr 


38 


68 


131 


s 


16 


24 


60 


Y 


39 


70 


134 


CI 


17 


26 


63 


Zr 


40 


72 


137 


Ar 


18 


27 


67 


Nb 


41 


74 


140 


K 


19 


29 


70 


Mo 


42 


77 


144 


Ca 


20 


30 


73 


Tc 


43 


79 


147 


Sc 


21 


32 


76 


Ru 


44 


81 


150 


Ti 


22 


34 


80 


Rh 


45 


83 


153 


V 


23 


36 


83 


Pd 


46 


86 


156 
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Table 2: The mass fractions of a-particles (X max ) and neutrons (X n ) out of which the seed 
nuclei are formed. 



Ye 


X max in % 


X n in % 


0.499 


99.8 


0.2 


0.498 


99.6 


0.4 


0.496 


99.2 


0.8 


0.49 


98 


2 


0.48 


96 


4 


0.47 


94 


6 


0.46 


92 


8 


0.45 


90 


10 


0.40 


80 


20 



Table 3: Abundances and mass fractions obtained for a boundary entropy beyond Sfi na i. 
The model parameters are S = 278, Y e = 0.49 and V exp = 7500 km/s (r exp = 35 ms). For 
further details, see text. 



Baryon 


Abundance (Y) 


Mass fraction (X in %) 


4 He 


0.2439 


97.55 


n 


0.22 lO" 1 


2.2 


H 


0.23 10~ 2 


0.23 


2 H 


0.28 lO" 6 


0.56 10~ 6 


3 H 


0.73 10- 7 


0.22 10" 6 


12 C 


0.22 10^ 8 


0.26 10~ 7 


16Q 


0.11 lO" 8 


0.18 10- 7 


3 He 


0.41 10~ 9 


0.12 10- 8 
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Table 4: The a-rich freeze-out conditions for an r-process at two expansion velocities V exp = 
7500 and 15000 km/s which correspond to expansion timescales of 35 and 15 ms, respectively. 
For a given Y e , Si n iu a i is the entropy which initiates an r-rpocess (Y n /Y see d ~ 1), Sfi na i is the 
last entropy which is able to build heavy seed nuclei (beyond it no r-process can take place), 
AS = Sfi na i — Siniuai, and ( ] represents the maximum Y n /Y seed ratio which can be 

V rseed / final 

attained when S = Sf ina i. For further details, see discussion in text. 



Y e 


7500 km/s 


15000 km/s 




^initial 


Sfinai 


AS | 


' Y n \ 

\ Yseed J final 


^initial 


Sfinai 


AS ( 


> Yseed J final 


0.499 


206 


253 


47 


17 


156 


193 


37 


17 


0.498 


185 


258 


73 


37 


143 


194 


51 


32 


0.496 


162 


263 


101 


60 


128 


198 


70 


50 


0.494 


153 


267 


114 


72 


121 


202 


81 


62 


0.492 


147 


272 


125 


84 


117 


205 


88 


68 


0.490 


142 


277 


135 


96 


113 


209 


96 


78 


0.486 


136 


283 


147 


109 


109 


216 


107 


92 


0.482 


131 


295 


164 


132 


105 


223 


118 


107 


0.478 


126 


303 


177 


148 


101 


229 


128 


120 


0.474 


122 


311 


189 


167 


98 


236 


138 


139 


0.470 


118 


319 


201 


188 


95 


242 


147 


155 


0.466 


114 


327 


213 


212 


92 


248 


156 


174 


0.462 


111 


334 


223 


233 


89 


254 


165 


195 


0.458 


107 


341 


234 


258 


87 


260 


173 


219 


0.454 


104 


348 


244 


286 


84 


265 


181 


239 


0.450 


100 


354 


254 


311 


81 


270 


189 


261 


0.400 


53 


416 


363 


755 


43 


320 


277 


655 
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Table 5: r-Process simulations with V exp = 7500 km/s and Y e = 0.45. The first column 
shows a choice of entropies, the second shows the corresponding r-process strength, i.e., the 
Y n /Y see d ratio. The third and sixth columns represent the r-process durations, i.e., the time 
between the freeze-out of the charged-particle reactions and the freeze-out of the following 
r-process. The fourth and seventh columns show the r-process freeze-out temperatures, and 
the fifth and last columns show the nucleus with the maximum abundance after decay back 
to stability for the FRDM and ETFSI-Q mass models, respectively. Note the earlier onset 
of fission cycling with FRDM. 



Entropy 


Strength 


FRDM 


ETFSI-Q 


S (ke/baryon) 




t freeze 


(ms) Tg 


Maximum 


t freeze 


(ms) T 9 


Maximum 


100 


1 


1 


2.95 


92 Zr 


1 


2.93 


92 Zr 


130 


6 


21 


2.14 


104 Ru 


34 


1.81 


102 Ru 


175 


23 


91 


1.09 


126 Te 


86 


1.13 


128 Te 


195 


34 


197 


0.626 


129 Xe 


138 


0.818 


130 Te 


215 


48 


327 


0.411 


132 Xe 


208 


0.600 


140 Ce 


230 


61 


375 


0.365 


132 Xe 


241 


0.531 


154 Sm 


236 


66 


398 


0.346 


132 Xe 


253 


0.511 


162 Dy 


260 


92 


528 


0.269 


195p t 


322 


0.417 


192 Os 


270 


107 


600 


0.239 


196p t 


398 


0.347 


195p t 


280 


121 


684 


0.212 


129 Xe 


470 


0.299 


206 Pb 


295 


146 


860 


0.171 


129 Xe 


586 


0.244 


232 Th 


300 


154 


926 


0.159 


129 Xe 


629 


0.229 


130 Te 
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Table 6: The main isobaric progenitors of 130 Te, 162 Dy and 195 Pt using the mass model 
ETFSI-Q. The model parameters are Y e = 0.45, V exp = 7500 km/s, and S = 195, 236 and 
270. The neutron, chemical and dynamical freeze-outs mean Y n /Y heavy < 1, the break-out of 
the (71,7) — (7, n)-equilibrium and \Y n (t) /Y n (t)\ > \p(t)/p(t)\, respectively. The superscripts 
b and e denote the "begin" and "end" of the (n, 7) — (7, n)-equilibrium, respectively. 



Peak 


Freeze-out 


Progenitor 


Y (10- 5 ) 


Time (ms) 


T 9 


n n (cm 3 ) 


Yn/Yheayy 


130 Te 


Neutron 


130p d 


50 


138 


0.82 


10 23 


1 




Chemical 6 


130 Pd 


23 


171 


0.71 


10 21 


10- 2 




Chemical 6 


130 Cd 


21 


414 


0.33 


10 18 


10- 4 




Dynamical 


130 Cd 


16 


415 


0.33 


10 18 


10- 4 


162 Dy 


Chemical 6 


162 Xe 


1.3 


220 


0.57 


10 24 


9 




Neutron 


162 Ba 


4.0 


252 


0.52 


10 23 


1 




Chemical 6 


162 Ce 


1.9 


331 


0.41 


10 19 


10- 4 




Dynamical 


i62 Nd 


0.7 


479 


0.29 


10 17 


10" 5 


195p t 


Chemical 6 


195 Hq 


1.9 10" 3 


212 


0.59 


10 24 


47 




Neutron 


195 Tm 


1.9 


396 


0.35 


10 22 


1 




Chemical 6 


195 Tm 


4.3 


422 


0.33 


10 21 


0.2 




Dynamical 


195 Yb 


12 


701 


0.21 


10 18 


10~ 4 



Table 7: The total mass of r-process material which can be ejected from the HEW of core- 
collapse supernovae for two different expansion velocities of V exp = 7500 and 15000 km/s. For 
each electron abundance Y e , the integration over the entropies starts from S initia i (Y n /Y seed rs 
1) to the maximum entropy Sfi na i(Y e ). For more details, see text. 





Mass in 10" 


" 4 M 


Y e 


7500 km/s 


15000 km/s 


0.450 


3.42 


3.30 


0.458 


2.71 


2.57 


0.462 


2.34 


2.28 


0.474 


1.43 


1.37 


0.482 


0.89 


0.85 


0.490 


0.43 


0.41 


0.494 


0.22 


0.21 


0.498 


0.05 


0.04 



High-Entropy Wind: V = 7500 km/s 
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Neutrons Seed 

a-Particles 




50 100 150 200 250 300 350 
Entropy, S(k B /baryon) 
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1e-05 
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HEW, FRDM, V exp = 7500 km/s, Y e = 0.45 
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